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I.  INTRODUCTION  AND  SUMMARY 


The  response  of  thin  targets  to  projectile  impact  is 
investigated  under  the  assumption  that  membrane  theory  is 
sufficient  to  describe  the  interaction.  Though  it  is  possible 
with  existing  computer  codes  to  compute  such  problems  in  , . 
greater  detail  as  described,  for  example,  by  Dienes  et  al’-1^, 
the  membrane  theory  is  amenable  to  analytic  treatment  since 
the  equation  of  motion  is  of  only  second  order.  This  makes 
it  possible  to  gain  an  overview  of  the  problem  which  exhibits 
the  main  dynamic  phenomena.  It  also  makes  it  possible  to 
analyze  approximately  the  response  of  materials  such  as 
cloth  and  skin  which  are  not  represented  by  bending  theory 
either  in  analytic  or  in  computer  solutions. 

The  organization  of  this  report  is  strongly  influenced 
by  a companion  document  by  Dienes  and  Miles (2)  which  describes 
the  membrane  theory  in  detail  in  a form  suitable  for  publi- 
cation in  the  open  literature.  In  this  report  the  main  con- 
clusions are  summarized,  and  applied  to  a variety  of  ballistic 
problems.  In  Section  II  the  one  dimensional  plastic  response 
of  a wire  to  impact  is  analyzed.  The  conclusions  are  similar 
to  those  for  axisymmetric  impact,  but  the  mathematics  is  con- 
siderably simpler  because  of  the  elementary  nature  of  solutions 
to  the  one  dimensional  wave  equation.  The  results  of  the 
axisymmetric  impact  problem  are  summarized  in  Section  III. 

One  of  the  main  conclusions  is  that  the  residual  velocity 
does  not  go  to  zero  as  the  impact  velocity  decreases  to 
the  ballistic  limit.  Projectiles  that  penetrate  will 
always  do  so  with  a substantial  velocity.  Residual  velocities 
of  penetrating  projectiles  are  compared  with  experimental 
data  in  Section  IV.  The  method  allows  data  for  a specific 
target  thickness,  projectile  size  and  projectile  velocity 
to  be  used  in  determining  the  ballistic  properties  under 
other  impact  conditions.  This  requires  that  a value  for 
the  ballistic  figure  of  merit,  w,  be  Obtained  for  each 
material.  Since  the  resistance  of  materials  to  impact  is 
determined  by  this  parameter,  it  cdn  be  used  to  rank  materials 
in  a quantitative  fashion,  even  when  the  tests  performed  are 
quite  dissimilar.  It  is  shown  in  Section  V that  this  approach 
works  successfully  for  a nylon  for  which  low  speed  data  is 
available.  At  high  speeds,  however,  penetration  at  the 
ballistic  limit  may  involve  cratering, .melting,  spallation 
and  other  phenomena  which  the  current  theory  does  not  allow 
for.  In  Section  VI,  the  figure  of  merit,  W,  is  given  for  a 
number  of  materials.  This  allows  them  to  be  ranked  in  order 
of  their  ballistic  figure  of  merit.  Finally,  conclusions 
and  recommendations  are  given  in  Section  VII. 
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In  Appendix  I it  is  shown  that  membrane  theory  cor- 
relates well  with  experimental  data  on  the  response  of  circular 
plates.  A listing  of  the  computer  program  used  to  evaluate 
the  ballistic  limit  and  residual  velocity  is  given  in  Appendix 
II. 


II.  ONE  DIMENSIONAL  IMPACT 


To  illustrate  the  physical  character  of  the  membrane 
approach  to  target  deformation,  consider  the  response  of  a 
rigid,  perfectly  plastic  wire  to  impact.  Since  the  tension 
in  the  wire  is  constant,  its  motion  is  governed  by  the  equation 

pAy  = aky"  (2.1) 


where  A denotes  the  cross-sectional  area  of  the  wire; 
y(x,t),  its  deflection;  p,  the  density  and  0 the  stress. 

In  view  of  the  assumption  of  ideal  plasticity 

a = a (2.2) 


where  Y is  the  flow  stress  and  z,  is  ±1,  the  stress  being 
positive  when  the  strain  rate  is  positive  and  negative  other- 
wise. Consequently,  the  equation  of  motion  reduces  to  the 
classical  wave  equation, 


• • 

y 


c2y ' 


(2.3) 


where 


c = /Y/p 


during  the  initial  phase  of  the  impact,  with  the  quantity,  c, 
being  interpreted  as  the  speed  of  plastic  waves. 

The  solution  can  be  written  in  the  form 


y = f(x  - ct)  + g(x  + ct) . 


(2.4) 


We  consider “only  the  right  going  waves  in  the  region  x > 0 
of  Fig.  1,  so  that  the  solution  reduces  to  the  form 


y = f (x  - ct)  . 


(2.5) 
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Fig.  1 Sketch  of  geometry  for  the  problem  of  an  infinite, 
ductile  wire  struck  by  a mass,  m. 


Only  the  region  to  the  right  of  the  projectile  need  be 
considered,  since  the  geometry  is  symmetric  with  respect  to 
the  center  of  the  impacting  mass.  The  boundary  condition 
to  be  imposed  at  the  origin  is 

mz  = hYy'  (2.6) 

where  m denotes  the  mass  per  unit  width  and 

z = y (o,t)  . (2.7) 

denotes  the  displacement  of  the  mass.  In  view  of  Eq.  1.5 
the  boundary  condition  can  be  written 

ay'  = y"  (2.8) 

where 


a = ha/c2m  = hp/m.  (2.9) 
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The  argument  of  y in  Eq.  1.8  is 


£ = - ct 


since  the  equation  is  valid  only  at  x = 0.  It  has  the 
solution 


f(5)  * y„(e“5  - 1) 


under  the  constraint  that  the  deflection  must  initially 
vanish.  It  is  straightforward  to  show  that  the  overall 
string  response  is,  then,  given  by 

y - y.Ce01'*  - ct)  - l)  . 


The  initial  condition  for  this  problem,  based  on  con- 
servation of  momentum,  is 


(m  +••  pha)y  = mv 


where  v is  the  initial  projectile  velocity.  Solving  for 
unknown  coefficient,  one  finds 


y = 


m 


m + pha 


m 

Kp 


1 - e 


a(x  - ct) 


) • x 


< ct 


= 0. 


X > ct 


At  late  times  the  deflection  approaches  the  value 


y 


00 


m m 
m + hpa  hpc 


v. 


Using  the  small  strain  approximation 


(2.10) 

(2.11) 

(2.12) 

(2.13) 

the 

(2.14a) 

(2.14b) 

(2.15) 

(2.16) 
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we  find  the  maximum  strain  equal  to 


- = 1 / m v \z 

max  1 \ m + pha  c ) 


If  it  is  assumed  that  the  material  is  ductile  and  fails  when 
the  strain  reaches  a critical  value,  e£,  (sometimes  termed 
the  breaking  index)  then  failure  occurs  when  the  impact 
velocity  exceeds  the  critical  value 


v = 


m + 


pha 


m 


✓2e£Y/p 


(2.18) 


The  quantity 


w = /2  e f Y/ p 


(2.19) 


is  characteristic  of  the  material  and  can  be  determined 
either  from  static  measurements  or  from  a measurement  of 
the  critical  velocity,  vc.  For  a rigid-plastic  material 
w2/2  is  equal  to  the  amount  of  energy  per  unit  mass  that 
the  material  can  absorb  without  failure. 
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III.  AXI SYMMETRIC  IMPACT 


The  response  of  an  infinite  target  impacted  by  a 
cylindrical  projectile,  under  the  assumption  that  the  flow 
stress  is  constant  in  the  target  material,  is  discussed 
in  detail  in  the  companion  report  by  Dienes  and  Miles.  (2) 
In  this  section  the  main  results  are  summarized  with  an 
emphasis  on  their  physical  interpretation. 

The  equation  of  motion  for  the  deflection > y,  of  a 
membrane  under  the  assumption  of  axial  symmetry  is 


pytt  = Y(ryr)r  (3.1) 

where  Y is  the  flow  stress;  p,  the  density  and  r,  the  radial 
coordinate.  On  the  contact  circle  the  membrane  stresses 
cause  a deceleration  of  the  projectile  which  leads  to  the 
boundary  condition 


2naYhy  = (m  + ira2ph)y  (3.2) 

r t u 

where  m denotes  the  projectile  mass;  a,  its  radius  and  h 
is  the  membrane  thickness.  It  proves  convenient  to  define 
the  dimensionless  variables 


n'  = r/a  (3.3) 

t = ct/a  (3.4) 


and 


? = y/(l  ' y)Ba 


(3.5) 


where 


y 


7ra2ph 
m + ita2ph 


(3.6) 
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and 


3 = v/c  .. 


(3.7) 


In  terms  of  these  dimensionless  variables  the  boundary 
value  problem  can  be  cast  into  the  form 


(r)£  ) = -n?  (3.8) 

*•  n n tt  ' 

with  the  boundary  condition  at  n = 1 


2«n  = (3-9> 

and  the  initial  conditions 


C = C = 0 (n  > 1,  x = 0) 

and 

= o,  i = i (t  = o)  . 


(3.10) 


(3.11) 


By  the  method  of  Laplace  transforms  it  can  be  shown 
that  the  solution  can  be  expressed  as 


t = 4y  Re  [AeS|>T|  + 2yJ(x) 


(3.12) 


where  Re  denotes  the  real  part  of  the  expression  in 
brackets.  Here 


A = 


1 

s3  + 4y(l  - y) s 
0 0 


(3.13) 


and 

i 
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J(T)  = 


/ 


e XTdx 


[x2K  (x)  + 2yK  (x)]2  + tt2[x2I  (x)  - 2yl  (x)  ] 2 (3.14) 

01  o i 


where  K?,  K,,  I0  and  Ia  are  modified  Bessel  functions.  In  the 
expression  for  A,  s0  denotes  the  root  of  the  transcendental 
equation 


sK  (s)  + 2yK  (s)  = 0,  (3.15) 

0 1 

$ 

The  term,  A,  arises  from  the  contribution  of  a pair  of 
poles  in  the  s plane.  The  maximum  dimensionless  deflection 
Smax>  was  obtained  numerically  and  is  shown  in  Fig.  2. 

i 

The  strain  in  the  membrane  is  given,  as  in  the 
previous  section,  by 


£ 


(3.16) 


for  small  deformations.  Denoting  by  a.  the  deceleration  of 
the  projectile,  it  can  be  readily  shown  that  the  strain  is 
given  by 


£ = \ (1  - u)  2 (3a) 2 • (3.17) 

Under  the  assumption  of  a critical  strain  failure  criterion, 
it  follows  that  penetration  occurs  for  impact  velocities 
above 


vc  = 


w 


(1  - U)  a, 


m 


(3.18) 


where  am  denotes  the  maximum  value  of  a.  As  in  the  previous 
section,  the  quantity  w can  be  taken  as  a figure  of  merit  and 
its  value  is  given  by 


w = flz  fY7p 


(3.19) 
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The  quantity  am  is  plotted  in  Fig.  3 as  a function  of  y. 

Its  value  always  exceeds  unity.  Thus,  the  ballistic  limit 
for  a cylinder  of  radius  a is  always  less  than  that  for  a 
rectangular  prism  of  width  a.  The  ballistic  limit  decreases 
as  y becomes  small  which,  from  inspection  of  (3.6),  occurs 
when  the  projectile  length  increases  and  the  radius  decreases 
with  a fixed  projectile  mass.  This  is,  of  course,  in  agree- 
ment with  one's  expectations.  The  dimensionless  time  at 
which  peak  deflection  occurs,  xd>  is  shown  in  Fig.  4,  which 
also  includes  the  dimensionless  time  to  maximum  strain,  xm. 
The  actual  times,  based  on  Eq.  3.4,  are 


t 


d 


atd/c 


(3.21) 


for  maximum  deflection,  and 


t 


m 


ax  /c 
m 


(3.22) 


for  maximum  strain.  Maximum  strain  always  occurs  prior  to 
maximum  deflection.  This  implies  that  when  failure  occurs, 
based  on  a maximum  strain  criterion,  that  the  projectile 
is  still  in  motion  and,  hence,  that  it  will  penetrate  at  a 
finite  velocity.  It  is  shown  in  Ref.  2 that  the  residual 
velocity  is 


vr  = F (y)  v. 


with  F(y)  being  plotted  in  Fig.  5. 


IV.  COMPARISON  OF  MEMBRANE  THEORY  WITH  EXPERIMENT 


It  is  known  that  for  small  deflections,  particularly 
when  the  materials  lie  in  the  elastic  regime,  the  deformation 
of  plates  is  governed  by  bending  theory.  With  increasing 
deflections,  stretching  of  the  middle  surface  becomes  important 
and  the  membrane  stresses  begin  to  dominate.  Ultimately,  for 
very  large  deflections,  the  resistance  is  governed  virtually 
in  its  entirety  by  the  membrane  stresses,  which  may  be  con- 
sidered uniform  through  the  thickness  and  independent  of 
radius.  A comparison  between  calculations  based  on  this 
approach  and  measured  deflections  of  impulsively  loaded 
circular  metal  plates  is  given  in  Appendix  I.  The  comparison 
shows  that  membrane  theory  describes  the  deformation  well  when 
the  deflection  exceeds  15%  of  the  radius. 

To  indicate  the  potential  accuracy  of  the  method  for 
determining  the  ballistic  limit,  we  consider  here  the  resistance 
of  a sheet  of  the  titanium  alloy  Ti5Al-2.5Sn.  The  residual 
velocities  measured  in  a series  of  tests  described  by 
Bruchey(3)  in  which  steel  cylinders  were  fired  into  sheets 
0.1295  cm  in  thickness  are  shown  in  Fig.  6 for  a cylinder  of 
radius  .284  cm  and  mass  1.037  gm.  Based  on  a titanium  density 
of  4.5  g/cm3  we  find  y = 0.125  and  F(y)  = 1.11.  If  we  take 
the  ballistic  limit  from  Fig.  6 as  the  lowest  striking 
velocity  at  which  residual  velocities  are  clustered,  about 
340  m/s,  then  w = 330  m/ s . Now,  the  flow  stress  for  a 
similar  alloy,  T16A1-4V,  is  given  by  Lindholm,  Yeakley  and 
Bessey  (4j  as  188ksi  (13.0  x- 109d/cm2)  and  the  cited  strain 
to  failure  based  on  post-test  measurements  is  18  percent. 

(The  similarity  in  properties  of  the  alloys  can  be  verified 
by  comparing  stress-strain  curves  given  by  Wolf (5)).  Using 
these  data  we  find  w = 324  m/s,  which  is  only  2 percent  less 
than  the  value  obtained  from  ballistic  data.  The  strain-to- 
failure  of  18%  based  on  post-test  measurements  is  signigicantly 
larger  than  is  indicated  by  their  stress-strain  curve,  but 
such  data  is  not  commonly  available  in  handbooks,  and  for  this 
reason  we  have  had  to  refer  to  the  more  specialized  data  for 
Ti6Al-4V. 

The  minimum  residual  velocity  can  be  obtained  using 
the  value  of  y cited  above  and  Fig.  7.  For  a projectile 
velocity  of  340  m/s  we  find  vr  = 210  m/s,  in  good  agreement 
with  the  lowest  values  observed  by  Bruchey.  This  result, 
that  the  residual  velocity  does  not  go  to  zero  as  the  impact 
velocity  decreases  to  the  ballistic  limit,  is  our  most 
striking  conclusion. 

To  determine  the  residual  velocity  above  the  ballistic 
limit,  we  note  that  at  r = a the  strain  is  given  by  Eq.  3.16. 
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The  value  of  yr  at  the  time  of  maximum  strain  may  be  obtained 
by  solving 


/2 IJ--  - (1  - u)0C(Tm)  (4.1) 

for  the  time  of  maximum  strain,  xm,  at  a specific  impact 
velocity  (recalling  that  < 0)  . Then  the  residual  velocity 
is  given  by 


vr  = 2y  (1  - u)vc(tm)  . (4.2) 

In  Fig.  7 we  have  plotted  ap  = -,C(xm)/2  versus  ?(2y)/2y 
for  various  values  of  y.  Thus,  s(xm)  can  be  obtained  from 
(4.1)  without  explicitly  determining  xm  by  means  of  the 
graphs.  In  practise,  it  is  convenient  to  determine  ap  from 
the  relation 


w/v 

af  2y (1  - y) 


(4.3) 


and  the  residual  velocity  is  given  by  4.2.  The  residual 
velocities  for  ap  _>  1 shown  in  Fig.  6 were  obtained  in  this 
way,  and  though  the  experimental  data  are  too  scattered  to 
confirm  the  theory  in  detail  in  this  range,  the  overall 
trend  is  reasonable.  As  the  impact  velocity  increases,  the 
time  to  maximum  strain  decreases  until  at  xm  = 0,  '£  = 1. 

The  corresponding  value  of  v,  which  we  denote  by  v,  is  the 
upper  bound  of  the.  impact  velocities  that  result  in  membrane 
deflection  followed  by  membrane  failure. 

At  impact  velocities  above  v the  membrane  fails 
instantly  (in  the  current  theoretical  model)  and  the  residual 
velocity  is  reduced  only  by  transfer  of  momentum  to  the 
ruptured  plug  that  is  removed  from  the  membrane,  Then 

vr  = !L—  v (4.3a) 

m + Tra2ph 


° U - y)v 


(4.3b) 


This  trend  agrees  well  with  the  measured  residual  velocities 
shown  in  Fig.  6. 
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V.  CORRELATION  OF  DATA  FOR  DIFFERENT  VELOCITIES 


The  theory  outlined  in  the  preceding  sections  is  based 
on  a number  of  strong  assumptions  concerning  projectile  shape 
and  rigidity  and  the  characteristics  of  the  target.  In  this 
section  we  consider  two  experimental  determinations  of  the 
ballistic  limit  and  examine  the  correlation  provided  by 
membrane  theory. 

Laible,  Figucia  and  Ferguson have  examined  the 
resistance  to  penetration  of  a number  of  high-modulus  fibers, 
primarily  nylons  and  other  organic  fibers.  In  one  pair  of 
experiments  the  ballistic  limit  (V50)  was  determined  for  an 
organic  fiber  described  as  X-500,  type  I.  Its  tenacity  is 
reported  as  12-14  grams  per  denier;  the  modulus  as  500  grams 
per  denier;  elongation  to  failure  as  2-4%;  and  density  as 
1.47  g/cm3.  The  relevance  of  these  properties  to  the 
properties  of  the  fabric  is  somewhat  obscure,  however, 
since  the  mechanical  behavior  of  a fabric  may  be  very  complex 
even  when  its  constituent  fibers  behave  in  a linear  fashion. 
The  tests  involved  impact  by  a 17  grain,  0.22  caliber  missile. 
Its  velocity  was  varied  until  the  value  at  which  50%  of 
the  projectiles  were  defeated  by  the  target  was  determined, 
and  this  value  was  defined  as  the  ballistic  limit,  V50.  The 
results  are  summarized  in  Table  I.  It  is  remarkable  that  the 
values  of  the  ballistic  figure  of  merit,  w,  agree  within 
2 percent  between  the  two  tests  in  view  of  the  many  respects 
in  which  the  details  of  the  encounter  differ  from  the 
assumptions  made  in  the  analysis. 

Bruchey(3)  has  determined  the  ballistic  limit  for  two 
steel  projectiles  fired  into  titanium  targets  of  the  same 
thickness.  The  results  for  the  heavier  projectile,  with  a 
mass  of  1.037  g,were  described  in  the  previous  section.  The 
projectile  radius  was  0.284  cm,  and  the  areal  density  of  the 
target  was  .583  g/cm2.  As  previously  described,  this  leads 
to  a value  of  330  m/sec  for  the  ballistic  figure  of  merit,  w. 

A cylinder  having  dimensions  which  were  approximately  half 
those  of  the  low  speed  projectile  (a  = .131  cm)  was  also 
tested,  and  the  ballistic  limit  was  determined  to  be  585  m/s. 
It  follows  that  the  value  of  y is  .1945,  F(y)  = 1.02  and 
w = 481m/s.  This  figure  is  significantly  higher  (46%)  than 
that  determined  from  the  low  velocity  projectile,  indicating 
that  at  higher  speeds  the  target  material  behaves  more 
efficiently.  This  can  probably  be  attributed  to  the  greater 
dissipation  of  energy  per  unit  mass  involved  with  high  speed 
projectiles,  with  a substantial  fraction  of  the  energy  going 
into  the  cratering  process.  It  would  be  of  interest  to 
examine  the  interaction  in  careful  experiments  or  numerical 
calculations  to  determine  in  greater  detail  why  the  target 
is  more  efficient  in  resisting  high  speed  impacts. 
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VI.  RANKING  OF  MATERIALS 


The  impact  resistance  of  a number  of  materials  for 
which  ballistic  data  are  available  is  shown  in  Table  II. 

Based  on  this  data,  the  Hadfield  steel  which  has  been 
traditionally  used  for  infantry  helmets  and  many  kinds  of 
armor  plate  ranks  lowest  among  the  materials  considered 
with  w = 260  m/s.  Kevlar,  recently  developed  by  Dupont, 
ranks  highest  among  the  materials  analyzed,  with  w = 872  m/s. 
These  data  are  taken  from  widely  differing  sources,  however, 
and  it  would  be  desirable  to  conduct  a systematic  series  of 
tests  to  provide  a more  direct  verification  of  this  approach 
to  ranking  materials.  It  is  not  entirely  clear  from  the  raw 
data  what  properties  constitute  a good  ballistic  material. 

It  may  be  that  woven  materials  are  better  than  plate  materials, 
and  it  may  be  that  local  melting  plays  a role  in  some  cases. 
Though  the  approach  proposed  here  may  serve  as  a rough  guide 
to  the  relative  efficiency  of  different  materials,  the  choice 
of  materials  for  specific  protective  missions  should  involve 
careful  consideration  of  the  specific  penetration  mechanisms. 
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COMPARISON  OF  THE  BALLISTIC  RESISTANCE  OF  VARIOUS 
MATERIALS  RANKED  BY  THE  PARAMETER  w 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  proved  possible  to  determine  the  motion  of  a 
membrane  impacted  by  a cylindrical  mass  in  closed  form  with  the 
assumption  that  the  flow  stress  in  the  membrane  material  is 
uniform.  A comparison  of  the  computed  deflection  with  measure- 
ments published  by  Florence  in  1966  has  indicated  that  the 
assumption  that  the  flow  stress  is  constant  and  equal  to  the 
yield  strength  of  the  material,  for  circular  plates, 
leads  to  reasonably  accurate  results  when  the  deflection 
exceeds  151  of  the  radius. 

Of  the  various  possible  failure  theories,  a critical 
strain  criterion  is  the  simplest  to  implement  in  the  current 
theoretical  approach.  It  appears  to  lead  to  good  correlation 
when  compared  with  the  low  speed  titanium  results  obtained 
by  Bruchey,  both  with  regard  to  the  ballistic  limit  and  the 
residual  velocity.  It  is  found  that  the  residual  velocity  is 
finite,  even  when  the  target  is  penetrated  at  the  ballistic 
limit.  This  is  due  to  the  fact  that  maximum  strain  occurs 
before  maximum  deflection  and,  consequently  that  failure 
occurs  before  the  motion  has  stopped.  In  the  example  of  the 
low -speed  titanium  impact  experiments  carried  out  by  Bruchey, 
the  theoretical  residual  velocity  was  210  m/s  at  the  ballistic 
limit  of  340  m/s. 

Configurations  for  which  the  stress  and  strain  are 
uniform  through  the  thickness,  as  assumed  in  membrane  theory, 
are  efficient  in  comparison  with  structures  that  undergo  bending 
stresses,  for  when  loads  are  resisted  by  bending  moments,  the 
middle  surface  is  strained  less  than  adjacent  surfaces, 
and  consequently  the  material  is  not  being  used  in  an  optimum 
fashion.  This  may  serve  as  an  explanation  of  a rather  sur- 
prising result  obtained  by  Laible  et  al. They  find  that 
the  addition  of  resin  to  fabric  laminates  lowers  the  ballistic 
limit.  In  one  example,  the  addition  of  8.5  oz/ft2  of  resin  to 
10.7  oz/ft2  of  fabric  lowered  the  ballistic  limit  from  830  to 
623  ft/sec.  This  suggests  that  the  filled  laminate  is  subject 
to  bending  stress  which  are  not  present  in  the  unfilled  laminate, 
and  that  failure  occurs  when  the  bending  stress  exceeds  the 
allowable  value.  This  occurs  at  a lower  velocity  than  when  the 
laminates  are  unbonded,  and  they  are  all  stressed  to  the  same 
level.  This  line  of  thinking  suggests  that  fabrics  are  more 
efficient  than  plates  of  the  same  material.  In  fact,  some 
tests  were  run  at  S3  in  which  woven  wire  (340  stainless)  was 
compared  with  plates  of  the  same  material.  The  woven  fabric 
turned  out  to  have  a relatively  low  ballistic  limit.  This 
was  subsequently  attributed  to  the  fact  that  the  wire  had 
an  elongation  of  only  151,  where  the  plate  had  an  elongation 
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of  401.  We  believe  that  if  a fabric  of  wire  which  was  not 
strain-hardened  were  fabricated,  then  it  would  indeed  prove 
more  efficient  than  plate. 

Few  materials  exhibit  a stress-strain  relation  which 
closely  resembles  the  rigid-plastic  behavior  assumed  in  the 
current  theory.  It  would  be  desirable  to  find  a good  method 
for  approximating  real  material  behavior  by  a solution  to 
the  linear  wave  equation.  That  is,  how  can  one  select  an 
equivalent  strength,  Ye,  which  in  some  sense  causes  the 
solutions  to 


py  = Y(y')V2y 


to  be  closely  approximated  by 
py  = Yev2y 

Procedures  for  this  kind  of  analysis,  termed  the  method  of 
equivalent  linearization,  have  been  developed  for  ordinary 
differential  equations,  as  discussed,  for  example,  by 
Booten(H)  an(j  Dienes, ^12),  but  they  do  not  appear  to  have 
been  developed  adequately  for  partial  differential  equations. 

The  problem  of  oblique  impact  has  not  been  solved,  even 
in  the  simple  case  of  a linear  membrane.  The  separation  of 
the  target  into  regions  of  tension  and  compression  seems  to 
present  a major  difficulty.  Perhaps  it  would  be  best  to 
use  a numerical  approach  to  some  of  these  problems,  retaining 
the  assumptions  of  membrane  theory,  in  order  to  gain  a better 
insight  into  the  mechanics  of  oblique  impact.  The  results 
might  suggest  some  approximations  that  would  be  useful  in 
the  oblique  problem. 
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APPENDIX  I.  RESPONSE  OF  CIRCULAR  PLATE  TO  UNIFORM  IMPULSE 


Florence (10)  has  determined  the  response  of  both 
aluminum  (6061-T6)  and  steel  (CR-1018)  circular  plates  to 
sheet  explosive,  which  provides  an  approximately  uniform 
impulse,  and  compared  the  maximum  deflection  with  that  cal- 
culated from  plate  (bending)  theory.  His  results,  which 
are  shown  in  Fig.  8,  suggest  that  bending  theory  is  inadequate 
for  deflections  that  exceed  101  of  the  plate  radius,  and  he 
remarked  that  membrane  stresses  might  account  for  the  dis- 
crepancies. We  proceed  to  apply  the  membrane  model  to  his 
configuration. 

The  notation  remains  as  in  Section  III,  except  that  a 
is  now  the  radius  of  the  circular  membrane.  The  motion  is 
governed  by  the  axisymmetric  wave  equation  (3.1),  the  boundary 
condition  (3.2)  is  replaced  by 


y = 0 (r  = a,  t > 0)  , (A.l) 


and  the  initial  conditions  are 


y ■ 0,  yt  **  vQ  s I/ph  (r  < a,  t = 0)  , (A.2a,b) 

where  I is  the  impulse  per  unit  area.  [In  reality,  vQ  is 
only  approximately  uniform  and  falls  sharply  to  zero  as 
r + a in  consequence  of  the  restraint  (A.l)  at  r = a;  however, 
this  deficiency  of  the  model  is  relatively  unimportant  for  the 
prediction  of  y in  r < a], 

• The  solution  of  (3.1),  (A.l)  and  (A. 2)  is  given  by 


y = ^ AiJ0(Xih)sin(XiT)  > 

i=l 


(A.  3) 


where  X±  is  the  i ' th  root  of  the  Bessel  function  Jq> 


Jn(X.)  = 0 (0  < X < X < ...<»)  , 

U 1 12 


(A.  4) 
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n and  t are  defined  by  (-3 -3)  and  (3.4),  and 


/v  J (A . n) ndn  o 

0 0 v 1 ' 2v  a 

A.  = . _o - o 

1 cAi  I*1  JQ2(Ain)ridn  cXi2Jj(Xi) 


(A.5) 


The  deflection  at  the  center  of  the  plate,  r = 0,  is  given 
by 


y0/a  = (vo/c)f(T)  , (A.  6>) 

00 

where  f(t)  = 2 S (A. 2J  (A . ) ] " 1 sin (A . t)  . (A. 7) 

Z- 4 1 l 1 1 

i=l 

Its  first  maximum  (determined  numerically)  of  1.0  is  the 
only  one  that  has  physical  significance,  since  the  formu- 
lation is  valid  only  for  increasing  deflection. 

The  theoretical  prediction  ymax/a  = l*°(vo/c) 
compared  with  Florence’s  results  in  Fig.  8.  It  is  evident 
that  the  membrane  model  provides  a significantly  better 
correlation  with  the  measured  results  than  does  the 
bending  theory  if  v0/c  > 0.15  and  that  the  correlation 
increases  with  v„/c  within  the  range  of  the  measurements. 
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APPENDIX  II 

LISTING  OF  COMPUTER  PROGRAM 
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END  OF  COMPILATION:-  NO  DIAGNOSTICS. 
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